storymap-info/process files.R

In Metropolitan-Council/planting.shade: Tree planting and preservation prioritization tool

library(sf)
library(tidyverse)
library(tigris)

######
# export for GEE
#####
nhood_list %>%  
  sf::st_write(., "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/neighborhoods.shp", append = FALSE)

#######
# exporting things for the storymap
######

#frogtown shapefile
nhood_list %>%
  filter(GEO_NAME == "Thomas-Dale/Frogtown") %>%
  sf::st_write(., "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/frogtown.shp", append = FALSE)

#st. paul neighborhoods
nhood_list %>%  
  filter(city == "St. Paul") %>%
  rename(Neighborhood = GEO_NAME) %>%
  sf::st_write(., "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/stpaul_nhoods.shp", append = FALSE)

#trees in tracts
mn_tracts %>%
  rename(tract_string = GEOID) %>%
  right_join(eva_data_main %>%
  filter(variable == "canopy_percent") %>%
  dplyr::select(raw_value, tract_string) %>%
  mutate(`Canopy coverage` = raw_value * 100)) %>%
  dplyr::select(tract_string, `Canopy coverage`) %>%
  rename(`Tract id` = tract_string) %>%
  sf::st_write(., "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/tree_tracts.shp", append = FALSE)

#bipoc in tracts
mn_tracts %>%
  rename(tract_string = GEOID) %>%
  right_join(eva_data_main %>%
               filter(variable == "pbipoc") %>%
               dplyr::select(raw_value, tract_string) %>%
               mutate(`Percent BIPOC` = raw_value * 100)) %>%
  dplyr::select(tract_string, `Percent BIPOC`) %>%
  rename(`Tract id` = tract_string) %>%
  sf::st_write(., "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/bipoc_tracts.shp", append = FALSE)


# #trees in blocks
# #blocks are just SO messy, it's not really worth it. Plus the census is unreliable at the block level anyhow
# #for some reason not working to dl off site: https://resources.gisdata.mn.gov/pub/gdrs/data/pub/us_mn_state_metc/society_census2010realign/
# blocks <- sf::read_sf(("~/Documents/GitHub/planting.shade/storymap-info/shapefiles/realign blocks/Census2010RealignBlock.shp")) %>%
#   filter(ALAND10 > 0,
#          Acres > 1) %>%
#   mutate(ratio = ALAND10 / AWATER10) %>%
#   filter(ratio > .05)
# blocks %>% arrange(ALAND10) %>% head() %>% data.frame()
# 
# leaflet() %>%
#   addTiles() %>%
#   addPolygons(data = filter(blocks, BLK10 == "270370608061012") %>% st_transform(4326))

# tree in block groups
tigris::block_groups(state = "MN", county = c("Ramsey", "Hennepin", "Dakota", "Carver", "Scott", "Anoka", "Washington"), year = 2010) %>%
  right_join(read_csv("./data-raw/TreeAcres/TreeAcres_blockgroups_year2020.csv", 
                      col_types = list("GEOID10" = "c"))) %>%
  transmute(tract_string = GEOID10,
            canopy_percent = `1` / ALAND10) %>%
  sf::st_write(., "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/trees_blockgroups.shp", append = FALSE)
  


usethis::use_data(ctu_list, overwrite = TRUE)




mn_tracts %>%
  rename(tract_string = GEOID) %>%
  right_join(eva_data_main %>%
               filter(variable == "canopy_percent") %>%
               dplyr::select(raw_value, tract_string) %>%
               mutate(`Canopy coverage` = raw_value * 100)) %>%
  dplyr::select(tract_string, `Canopy coverage`) %>%
  rename(`Tract id` = tract_string) %>%
  



###########
#export native lands
###########
oceti <- jsonlite::fromJSON("https://native-land.ca/api/index.php?maps=territories&name=oceti-sakowin-sioux,/indigenousTerritories.json") 
st_as_sfc(oceti, GeoJSON = TRUE)
meuse_sf = st_as_sf(oceti, coords = c("x", "y"), crs = 28992, agr = "constant")
st_write(meuse, "meuse.geojson")


rgdal::writeOGR(oceti,dsn= "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/oceti.GeoJSON", layer="oceti", driver="GeoJSON")

# oceti_geoms = oceti$geometry
# oceti_types = oceti_geoms$type
# oceti_geoms = oceti_geoms$coordinates
# oceti_plgs = lapply(1:length(oceti_geoms), function(x) {
#   iter_type = glue::glue("{toupper(oceti_types[x])}")
#   iter_geom = oceti_geoms[[x]]
#   iter_geom = glue::glue("{iter_geom[,,1]} {iter_geom[,,2]}")
#   iter_geom = paste(iter_geom, collapse = ', ')
#   iter_out = as.character(glue::glue("{iter_type}(({iter_geom}))"))
#   iter_out
# })
# oceti_properties = oceti$properties
# oceti_properties$geometry = oceti_plgs
# oceti_properties = sf::st_as_sf(oceti_properties, wkt = "geometry")

jsonlite::write_json(oceti, "~/Documents/GitHub/planting.shade/storymap-info/shapefiles/oceti.json")

read_sf("https://native-land.ca/api/index.php?maps=territories&name=oceti-sakowin-sioux,/indigenousTerritories.json") 





########
#city tree canopy
######
st_drop_geometry(ctu_list) %>% 
  full_join(st_drop_geometry(nhood_list)) %>%
  dplyr::select(GEO_NAME, canopy_percent) %>%
  filter(GEO_NAME != "Minneapolis",
         GEO_NAME != "St. Paul") %>%
  arrange(canopy_percent)


st_drop_geometry(nhood_list) %>%
  filter(city == "St. Paul") %>%
  dplyr::select(GEO_NAME, canopy_percent) %>%
  arrange(-canopy_percent) %>%    # First sort by val. This sort the dataframe but NOT the factor levels
  mutate(GEO_NAME=factor(GEO_NAME, levels=GEO_NAME)) %>%   # This trick update the factor levels
  ggplot( aes(x=GEO_NAME, y=canopy_percent)) +
  geom_segment( aes(xend=GEO_NAME, yend=0)) +
  geom_point( size=4, color=councilR::colors$cdGreen) +
  coord_flip() +
  theme_minimal() +
  xlab("")


  mutate(GEO_NAME = fct_reorder(GEO_NAME, -canopy_percent)) %>%
  ggplot(aes(x = canopy_percent,
             y = GEO_NAME)) +
  geom_segment(aes(yend=canopy_percent, yend=0)) +
  geom_point(size = 4)